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The finite size effects of the nanoparticles to the atomic pair distribution func- 
tions (PDF) are discussed by calculating the radial distribution functions (RDF) 
on nanoparticles with various shapes, such as sheet, belt, rod, tube and sphere, 
assuming continua. Their characteristics are shown depending on the shapes and 
the sizes of the nanoparticles. Alternately, these PDFs can be used to measure 
the shapes and the sizes of ordered lattice part inside of any materials such as 
nanoparticles and bulks. 
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1. Introduction 

Many types of nanoparticles have been synthesized, and stud- 
ied on their physical and chemical properties. The nanoparticles 
may exhibit the exotic physical properties, which can be differ- 
ent from those of the bulk materials due to their finite sizes. 
The studies on the atomic-scale structures of the nanoparticles 
give us information to understand their properties. However, it 
is hard to determine their structures by ordinary techniques used 
for the bulk materials, for example, traditional Rietveld anal- 
ysis. In the case of the nanoparticles, the well-defined Bragg 
peaks are not observed and the diffraction data mainly consist 
of the diffuse scattering, because the periodicity of their unit 
cells is limited to nanoscale. 

The technique of the atomic pair distribution function (PDF) 
can be applicable to the determination of the local structure of 
the nanoparticles. For the system composed of one kind of atom, 
the number of pair-atoms in a shell of thickness dr at distance 
r from another one are obtained as R(r)dr, where R(r) is the 
radial distribution function (RDF). R(r) is related to the reduced 
pair distribution function G(r) via the pair distribution function 
g(r) as follows (Egami & Billinge, 2003). 



g(r) = R(r)/4Trr 2 p , 
G(r) = 4wrpo[g(r) - 1] 



(1) 

(2) 



where po is a number density of atoms in the sample. Experi- 
mentally, G(r) can be obtained from the total scattering struc- 
ture function S(Q) via Fourier transformation as follows. 



G(r) 



Q[S{Q) - l}sin(Qr)dQ, 



(3) 



where Q is the magnitude of the wave vector. In the cases of 
amorphous materials and the nanoparticles, even if S(Q) does 
not have any well-difined peak as mentioned above, G(r) can 
have sharp peaks, at least, in the small r-region. Then the PDF 
analysis is applicable to the structural analysis of the materi- 
als whose structures have short range correlation. Recently, the 
technique of PDF has been applied for the determination of the 
structure of nanoparticles (Gilbert et al., 2004, Petkov et ai, 
2004, Gateshki et al, 2004, McKenzie et al, 1992). 



So far, the effect of the shape and finite size of a nanoparti- 
cle, however, has not been considered in the PDF analysis. The 
atomic correlation remains only in the size, and as a result, G(r) 
does not have any peak in the larger r-region than the particle- 
size if there is no correlation between near-by particles. Actu- 
ally, G{r) of Ceo has sharp peaks at r smaller than the diam- 
eter of C6o molecule while at r larger than the diameter, G(r) 
has only small broad peaks corresponding with the correlation 
between the molecules (Egami & Billinge, 2003). Furthermore, 
the intensity of the PDF must be reduced from that of the bulk 
sample with infinite size, due to the finite size even at r smaller 
than the particle size. So the distribution functions modified by 
their sizes and shapes need to be used for the detailed structural 
analyses of nanoparticles. 

Such analysis gives another structural information of the 
nanoparticle. The nanoparticle of the zinc sulfide is analyzed by 
PDF with the spherical shape effect (Gilbert et al, 2004). After 
they determined the averaged size and the shape of the nanopar- 
ticle by means of small angle X-ray scattering and ultraviolet- 
visible absorption spectroscopy, they discussed the disorder and 
the strain in the nanoparticle from the difference between the 
averaged size and the local correlation size by PDF. The correc- 
tions of the averaged size are crucial to discuss the local lattice 
disorder in the nanoparticle. 

In this paper, we calculate the RDFs of the various nanopar- 
ticles by assuming that they are continua, for the corrections of 
the averaged size and the shape effect, and the correction fac- 
tor to the RDFs and reduced PDFs for the various nanoparticles 
are obtained. By using the correction factor, the corrected for- 
mulation of the PDF analysis on nanoparticles are presented. 
Based on our calculations, the method to estimate the sizes of 
the parts with ordered and disordered lattice in the particle, is 
also proposed. 

2. Calculations 

In the calculation of the radial distribution function R(r) of a 
nanoparticle, we consider the atomic pair distribution only in 
the particle, and assume that the atomic density is zero (vac- 
uum) outside of the particle. The total atomic density in the 
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nanoparticle is p' Q . The RDF of the three dimensional continuum 
with infinite size is given as Roo (r) = A-Kr 2 p' Q , as mentioned in 
the next section. So the RDF of the nanoparticle, R nano (r), is 
modified by the correction factor f(r) which is defined as, 



fir) =R nano {r)IR OQ {r). 



(4) 



When r — > 0, f{r) must be unity, and when r is larger than the 
size of the nanoparticle or r — > oo, f(r) becomes zero. This 
factor can be regarded as a kind of particle form factor instead 
of a well-known atomic form factor. 

In the real analysis, the pair distribution between the particles 
must be considered. The formulation including the pair distri- 
bution beween particles will be discussed in §4. 

3. Calculated Results 

3.1. Simple Case 

In this subsection, we take up simple cases and discuss the 
effect of the dimensionality of particles to the radial distribu- 
tion function. First, we consider the wire with an infinitesimal 
thickness shown in Fig. 1(a). 



(c) 










ML 





Figure 1 

Schematic diagram of (a) ID, (b) 2D, (c) 3D particles, and (d) a nanobelt. 

The density in this wire is defined as p' Q lD whose dimen- 
sion is an inverse of the length. Here we assume that r is much 
smaller than the length of the wire L. The number of atoms in 
an length dr at the distance r from the atom at origin (O), is 
2p' Q 1D dr, as shown in the figure. Then the radial distribution 
function is obtained as, 



R w (r) = 2p{ Y 



ID- 



(5) 



In the case of the sheet with an infinitesimal thickness shown in 
Fig. 1(b), if r « L, atoms paired with the one at origin O, are 
in the ring with a radius r and a thickness dr which is shown by 
the bold line in Fig. 1(b). So Rwir) is given by, 



RiD{r) = 2irrp' Qt2D . 



(6) 



In the case of the three dimensional block shown in Fig. 1(c), 
paired atoms are in the spherical shell with a radius r and the 
thickness dr, and /?3o(r) is given by, 



R3D{r) = 47rr 2 /?o.3D- 



(7) 



In the cases of one, two and three dimensional nanoparticles, 
RDF functions have the r-dependences of r-constant, r-linear 
and r-square, respectively. 

Next, we calculate the RDF of the nanobelt with the width 
of a and an infinitesimal thickness which is shown in Fig. 1(d), 
as an example of the typical nanoparticle. Here, it is assumed 
that r <C L. In this case, we consider the number of atoms in 
an annulus of thickness dr at a distance r from another one at 
a position P which is distant from the edge of the belt by p, as 
shown in the figure. It is given by R(r 1 p)dr. So R(r,p) can be 
regarded as the "partial radial distribution function" (PRDF) at 
p. R(r) is given by J ( " R(r, p)dpla. If r < a/2, R(r, p) is propor- 
tional to the length of the circumference for r < p < a — r, and 
it is proportional to the length of the arcs with interior angles of 
2ir — 29 for p < r, and 2ir — 29' for a — r < p < a, where 
9 = cos" 1 £ and 9' = cos" 1 So 



RlDbeltir) 



(2ir - 29)rdp 



[ (2ir - 29')rdp 

J a — r 



2irrdp 



P0,2D la - 



In the case that a/2 < r < a, R(r, p) is proportional to the length 
of the arcs with interior angles of 2tt — 29 for p < a — r and 
2tt - 29' for r < p < a. At the other p, R(r,p) is proportional 
to the length of the bold line shown in Fig. 1(d). Then 



RlDbeltir) 



[2ir-29)rdp+ / (2tt - 29 - 29')rdp 



+ (2tt- 29')rdp 



Po,2D /a - 



Since for r > a, R(r, p) is obtained by the bold line for all p, 



RlDbeltir) = 



[ (2n - 26 - 26')rdp 
Jo 



P0,lD /a - 



By calculating the above integrations, we obtain the RDF of the 
nanobelt as follows. 



RlDbeltir) = i2nar - Ar 2 )p^ 2D la (r < a), 



R 



IDbell 



2 -1 a 

2-rrar — 4r — 4arcos - 



+4r\/ r 2 — a 2 p' 02 D^ a 



(8) 



When r is much smaller than a, RiDbeit ~ 2nrp' 2D , which is 
equal to Rmir). For r » a, RiDbeit ~ 2ap' 2D , equal to Rwir). 
The calculated /?(r) of the nanobelt is shown in Fig. 2. 
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Figure 2 

The radial distribution functions of the wire (ID), the sheet (2D) and the 
nanobelt. 

Here, we use parameters of p' 02D = 1. The RDFs of the 
wire and the sheet are also plotted in the figure, assuming same 
atomic densities. For r — > 0, the RDF of the nanobelt approach 
the RDF of the 2D sheet, and for r » a, they merge into that of 
the wire. We know the effect of the dimensionality to the RDF 
by such simple calculations. 

3.2. Nanosheet 

From this subsection, we consider the realistic model. First, 
we take up the nanosheet with a thickness of t as shown in Fig. 
3. 



9 = cos" 1 f and 9' = cos" 1 

Rsheet (r) 



. Then, for r < t/2, 

■t—r 



I 



j 2nr 2 (l + ^ dp + J 4-nr 2 dp 



2lT r 2 { 1 + *—E ] dp 



In the case of t/2 < r < t, R(r,p) = 2nr 2 (l + p/r)p' 03D for 
p < t - r, and R(r,p) = 2nr 2 {\ + '^p'^o for r < p < t. 
For t — r < p < r, R{r,p) is given as a product of p' 03D and 
the surface area of the object obtained by rotating the bold line 
shown in Fig. 3(b). For t/2 < r < t , 



2irrtdp 



Rshee,{r)= 2nr 2 (l + ^dp 

+ / 2 ^( , + £ 7 £ )* 



Since in the case of r > t , R(r, p) = 2irrt for all p, 

Rsheet {r) = / 2nrtdpp' 03D /t. 
Jo 

By calculating the above integrations, we get the RDF of the 
nanosheet as follows. 

R s heet(r) = (47rr 2 f - 2irr 3 )p' Q3D /t (r < t) 
Rsheet(r) = 2irtrp' QJD (r>t). (9) 



(a) 




(b) 




Figure 3 

(a) Three dimensional scheme of the nanosheet with a thickness t. (b) The sec- 
tioned diagram of the nanosheet. 

The area isLxL and it is assumed that r <C L. The atomic 
density in the sheet is p' Q 3D . Here, we also consider the PRDF at 
p, R(r, p), where p is a distance between the center of the sphere 
with a radius r and the bottom wall of the sheet, as shown in Fig. 
3(b). R(r,p) is obtained as product of p' Q3D and the volume of 
the overlapping part of the sheet and the spherical shell with a 
radius r and a thickness dr, shown by the bold line in the fig- 
ure. R s heet{r) is given by J' R{r, p)dplt. In the case of r < t/2, 
R(r, p) is given as the number of atoms in the complete spheri- 
cal shell for r < p < t — r. R{r, p) is proportional to the surface 
area of the object obtained by rotating the fan with a interior 
angle 9 and 9' for p < r and t — r < p < t, respectively, where 




Figure 4 

(a) The radial distribution function and (b) the correction factor of the nanosheet 
with thicknesses t. 

Figure 4(a) shows the RDF calculated on the nanosheet. In 
the calculations, the atomic density in the nanosheet p' 03D is 
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unity. For r -C t, the RDFs are proportional to r 2 . This r- 
dependence corresponds with the case of three dimensional par- 
ticle mentioned in §2. For r > t, R s heet{r) oc r, corresponding 
with Rio{r). The correction factor f sheet (r) which is defined by 
eq. (4) is shown in Fig. 4(b). At r = 0, f{r) is unity, and it 
linearly decreases with r in the region < r < t . 



are given by the relations, 



<j>i{x) = sin 
02 (x) = sin 



\Jr 2 — x 2 
-l a-p 
\J r 2 — x 2 



3.3. Nanobelt 

Next, we consider the case of the nanobelt with a width a and 
a thickness t , as shown in Fig. 5. 
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Figure 5 

(a) Three dimensional scheme of nanobelt with width a and thickness t. (b) The 
sectioned diagram of the nanobelt. (c) The diagram viewed along ^-direction. 

The length of the belt L is much larger than r. Figure 5(b) 
shows a sectioned diagram of the belt. The center of the cir- 
cle (P) in the figure is distant by p from the left wall and by q 
from the bottom wall. The PRDF at (p,q), R(r,p,q), is pro- 
portional to the surface area of the part of the sphere which 
cross the belt. The surface area of the object obtained by rotat- 
ing by <j)(x) around x axis with the radius y(x) is obtained by 
S = J y(x)d(j)ds, where ds = 'dx 2 + dy 2 . Since in the 

case of sphere, y = \J r 2 - x 2 , S = J r4>(x)dx. Here, we calcu- 
late R(r,p, q) by separating six parts as shown in Fig. 5(b). The 
surface area of L2 is obtained by the rotation from zero to <j>\ (x) 
and from ir — <p\ (x) to ir, and the area of R2 is obtained by the 
rotation from zero to 4>i{x) and from 7r — <pi(x) to n, as shown 
in Fig. 5(c). The surface areas of the other parts are obtained by 
the rotation of w. Then R(r,p, q) is given by 

R(r,p,q) =Po,3£> / nrdx+ I 2(f>\(x)rdx + I irrdx 
Uli JL2 Jls 

+ / nrdx + / 2(fi 2 (x)rdx+ / nrc, 

JR\ JR2 JR3 



(10) 



The direction of x-axis is shown in Fig. 5(b), and the origin of pr- 
axis corresponds with the center of the sphere. <j>\ (x) and fa (x) 



as shown in Fig. 5(c). The integral ranges in eq. (10) depend 
on the relations between r, t, a, p and q. The integral ranges of 
the left side (LI, L2 and L3) for various conditions of the above 
parameters are shown in Tabel 1 . 

Table 1 

Integral ranges of the left side (LI, L2 and L3) for various conditions of the 
relations between r, t, a, p and q. 



conditions 




LI 


L2 


L3 


(i) r < p 
r < q < t — r 
q < t — r, q > r 
q < t — r, q < r 
t — r < q < r 




r dx 

1-7 dx 

f-gdx 

[''"dx 

J -q 






(ii) r > p 



t - s/r 2 - p 2 < q < r 2 - 
\Jr 2 — p 2 < q < r and, 
q>t- v / r 2 -p 2 



f-9 dx 

J —Q 



' dx f 



q > r, q > t 



rf^dx r 



q < \/ r 2 — p 2 and, 
q < t - r 



V r 2 — p : 
-« 



. dx 



q < \J r 2 — p 2 and, 
t-q<q<t~y/r 2 - 

r < q < t — r 
q > r and, 

t — r<q<t — \/r 2 - 1 
\J r 2 — p 2 < q < r, and 
q < t - r 

^ ' r 2 — p 2 < q < r and, 
t-r<q<t-y , r 2 - l 



dx J 



si^dx 




— ~dx 

-p 1 


J Vr2 - p2 dx 


rt-q 


dx 

-P 2 

— ~dx 
-p 


r\A 2 -P 2 d 


rt-q 

W' 2 - 


-p 2 






— ^dx 

-p 1 




rt-q 

V'- 2 - 


dx 

-p 2 



For Rl, R2 and R3, the integral ranges are obtained by sub- 
stituting a — p for p in the table. The case shown in Fig. 5(b) 
corresponds with the condition of r > p, \J 1 r 2 — p 2 < q < r 
and t — r < q < t — \J r 2 — p 2 . Because the integrations for 
the range L2 and R2 can not be calculated analytically, they are 
obtained by the numerical calculations. Rb e it{ r ) is given by, 

Rbeit(r) = / / R{r,p,q)dqdpl(ta). 
Jo Jo 

The RDFs calculated for the nanobelts with various thickness t 
are shown in Fig. 6(a). 
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Figure 6 

(a) The radial distribution functions and (b) the correction factors of the 
nanobelts with various t. 

For r » ci,f, Rbeit(r) become flat and they are proportional 
to atp' 3D , corresponding with Rw(r). The correction factors 
obtained by eq. (4) are shown in Fig. 6(b). At r <~ 0, fbeit{r) 
for t = a and t = a/2 are slighty larger than unity, due to the 
insufficiency of the accuracy of the numerical calculation. 

3.4. Nanorod and Nanotube 

In this subsection, first, we consider the nanorod as shown in 
Fig. 7(a). 



(a) 






Figure 7 

(a) Nanorod with radius a. (b) The sectioned diagram of the nanorod whose 
center is located on A. (c) The diagram viewed along x-direction. 



We also assume that r <C L. Figure 7(b) shows a sectioned 
diagram. We consider the surface area of the overlapping part of 
the sphere with a radius r and the rod with a radius a. The cen- 
ter of the sphere is distant from the center of the rod by distance 
PA=p. The PRDF at p, R(r, p) is proportional to the surface 
area. It is similar to the cases of the nanosheet and the nanobelt. 
R(r, p) is obtained as 



/ 4(j>(x)rdx + / 2-rrrdx 
Jn hi 



{x-pf 



R{r,p) = 
where <f>(x) is given by 

4>{x) = sin - 



as shown in Fig. 7(c). The integral ranges II and 12 depend on 
the relation between the parameters, a, r and p. First, the case 
of r < a is considered. In this case, when p — a < —r, the 
sphere is perfectly enveloped in the rod. If p — a > —r, the 
integral ranges on x are obtained as, II:/? — a < x < r cos a and 



I2:rcosa < x < r, where cos a 
given for r < a, 

R(r, p) = 4i:r 2 p' Q 3D (p < a - r) 



p 2 +r 2 - 
2pr 



R(r,p) 



4<f>(x 



)rdx+ I 



2nrdx 



Then R(r,p) is 



P03D 



4r 



sin 



p—a 



(x - p) : 



dx 



+2nr 2 1 - 



,2\ I 



2pr 



P0.3D 



(p>a- r). 

(11) 



Since the RDF is given by R(r, p)2TTpdp/na 2 , R ro d{r) is 
obtained from the above equations as 



pa—r 

Rmd(r) = p' Q3D / 4irr 2 x lirpdp/- 
Jo 

i)3D / 

J a—r 



+ Po,: 



2irr z 1 



9 9" 

/r + r — a' 



2pr 



2irp 



+4r 



/ sin -1 

J p—a 



(x - p)< 



dx 



x 2-Kpdplira 2 (r <a). 



(12) 



The integration on x must be obtained by the numerical calcu- 
lation. 

In the case that r > a, R(r,p) can be given by the second 
formula in eq. (1 1) when p + a > r. When p + a < r, the circle 
with a radius a which corresponds with a section of the rod, is 
completely enveloped by the circle with a radius r. Then R(r, p) 
is given by 

rp+a 



rp+a 

R (r,P) = Po,2D / 4<P(x)rdx 

J p—a 
= 4r P03D 



p—a 



. _! a 2 - {x-p) 2 
sin \ ^ ^-t—dx. 



(13) 
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As the results, for r > a, R roc i is obtained as 



rr— a rp+a 

Rmd{r) = Po,3D / / 4rsin~ 

JO Jp-a 

x 2irpdp/na 2 



a 2 — (x — p) 2 



dx 



' r 

0,3D / 

J r—a 

f 



2,r 2 [l- pl + r2 - a2 



2pr 



4r I sin 

' p—a 

-2 



a 2 — (x — p) 2 



dx 



x 2tt pdp/ira (r > a). 



(14) 



The integrations on x included in the first and the third terms are 
calculated numerically. 

Next, we consider the nanotube shown in Fig. 8. 



(a) 






Figure 9 

(a) The radial distribution functions and (b) the correction factors of the nanorod 
and the nanotubes with various t. 



Figure 8 

(a) Nanotube with an external diameter 2a and an internal diameter 2b. (b) Sec- 
tioned diagram of the nanotube. 

An external diameter and an internal diameter are 2a and 2b, 
respectively. In the calculation, we consider the PRDFs at p for 
the nanorods with radii a and b, R a (r, p) and Rb(r, p), where p 
is a distance between the center of the tube and the sphere with 
a radius r. The RDF of the nanotube is given by using R a (r, p) 
and Rb{r, p), as follows. 

Rmbe(r)= f [R a (r,p)-R h (r,p)}2irpdp/ir(a 2 -b 2 ). (15) 

J b 

R a (r,p) can be calculated by same way as the case of the rod. 
Rb(r, p) can be given for the conditions of the parameters, 



R h (r,p) =4rp' 03D 

(p<r-b), 
R b (r,p) = p' 03D 

rn-b 



p+b 

p-b 



sin \ , = dx 



2^1 1-* 2 -^-"' 
2pr 



+Ar 



J 

J r c 



sin 



b 2 - (x - p) : 



dx 



R b (r,p) = ( P >r + b). 



(P>r-b), 
(16) 



By using the above equations, we calculate the RDFs of the 
nanorod and the nanotubes, as shown in Fig. 9(a). 



In these figures, a thickness t = a — b. The atomic densities 
in the rod and the tube, p' 3D are unity. At r ^> a,t, the RDFs 
have constant values which correspond with the a product of 
the sectioned area and the atomic density. It is same as the case 
of the nanobelt. R tu he{r) with a thin thickness have sharp peaks 
at r <~ 2a. The correction factors f tu be{r) are shown in Fig. 9(b). 
The correction factor of the tube decreases with r, and has a 
shoulder at r <~ 2a — t . 

3.5. Sphere 

In this subsection, the RDFs for the sphere and the spherical 
shell are calculated. Fisrt, the filled sphere with a radius a as 
shown in Fig. 10(a) is considered. 





Figure 10 

Sectioned diagram of (a) filled sphere and (b) spherical shell with external and 
internal radii a and b, respectively. 



research papers 



In this case, we also consider R(r,p) where p is a length PA 
shown in Fig. 10. In the case of r < a, it can be obtained as 

R(r, p) = 4nr 2 p' 03D {p<a-r) 



R ( r ,P) = Pq,3d I 



2tt rdx 



27T r 2 p' 3 D 



p 2 + r 2 — a 2 
2pr 



(p>a-r). 

(17) 



In the case of a < r < 2a, 
R(r,p) = (p<r-a) 
R ( r ,P) = P[)3D / 2lTrdx 

J r cos a 



p 2 + r 2 - a 2 



(p >r-a). 



For r > 2a, R(r, p) = 0. The RDF can be given by 

^R(r,p)4TT P 2 dp 



RFspherei,^} 



4 3 

^ira i 



(19) 



The above integrations of eqs. (17) and (18) give same results. 
Then 



RFsphere(r) = ^r 2 p' Q3D 



1 /r\3 
Rf sphere^) = (r>2a). 



/r\ 3 r 
- -3- 
\aJ a 



+ 4 



(r < 2a), 



(20) 



It is consistent with the result obtained by Mason. (1968) 

In the case of the spherical shell with a thickness t = a — b 
shown in Fig. 10(b), we consider the PRDF at p for the filled 
spheres with radii a and b, R a (r,p) and Rb(r,p), respectively, 
and the RDF can be given by 



REsphere( r ) 



J*lR a (r,p)-R b (r,p)}4Trp 2 dp 



f 7T (fl 3 - b 3 ) 



(21) 



R a (r,p) corresponds with eqs. (17) and (18). Rh(r, p) is obtained 

as 

R b {r,p)=0 (p<r-b) 

2 J 



R b {r,p)=2nr 1 p'^ D ^1 - 

(r-b<p<r + b) 
R b (r,p)=0 (p>r + b) 



2pr 



(22) 



By using eqs. (17), (18), (21) and (22), the RDFs for the spher- 
ical shells with various t are calculated by using ,00,30 = 1, as 
shown in Fig. 11(a). The correction factors are shown in Fig. 
11(c). 
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Figure 11 

(a) The radial distribution functions and (b) the correction factors of the 
nanospheres with various t. 

4. Discussions 

In this section, we present the formulation of the PDF analysis 
which takes account of the finite size effect of the nanoparti- 
cle. In the previous section, we neglect the atomic pair distribu- 
tion between the particles, however it must be considered in the 
real analysis. Here, we consider the partial raidal distribution 
function including the atomic pair distributions inside of a par- 
ticle and between the particles. It is assumed that the measured 
nanoparticles are in a vacuum (there is not any scatterer or atom 
outside of particles). Because in the real case, each particle is 
usually distributed randomly and there may be no atomic corre- 
lation between particles, the atomic distribution outside of the 
particle can be regarded as being continuous with a density of 
po which corresponds with the density averaged in the whole 
sample including all particles. The atomic arrangement inside 
of particles is periodic with a density of p' Q . Here we rewrite 
the radial distribution and the partial radial distribution func- 
tions obtained for the single particle made of atoms as R„ an o(r) 
and R„ano{r, p) (or R nan o{r, p, q)), respectively. Then the partial 
radial distribution function R(r,p) (or R(r,p,q)) are given as 
follows. 

R{r,p) = R na no{r,p) + [4tit 2 - S„ a „ (r,p)] po 

where S nano {r,p) shows a surface area that the sphere with a 
radius r centered at P has inside of the particle. In the above 
equation, the first term shows the pair distribution inside of 
the particle and the second term shows the pair distribution of 
outside of the particle (the pair distribution between the par- 
ticles). If the nanoparticle consists of continua, R„ a „ (r, p)l 'p' 
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becomes equal to S nano (r,p). Then the radial distribution func- 
tion is given by 

R(r) = R„ano{r) + [47rr 2 - S nano (r)] p Q , 

where S nano {r) is the surface area inside of the particle averaged 
by p. In the previous section, we know that 47rr 2 /(r) is substi- 
tuted for S nano (r). Then the RDF can be represented as 

R(r) = /(rjRoo (r) + 47rr 2 P o [1 - /(r)] , (23) 

where /?oo (r) is RDF of the bulk sample with infinite size. Here 
we use the relation 

f(r) = Rnano{r)IRoo{r), 

for the first term. For the bulk sample, the PDF and the reduced 
PDF are given by 

gocO) = R oo {r)/4Tvr 2 p' 
G oo ( r )=4wrp' [g oo (r)-l}, (24) 

where goo (r) and G^ (r) are the PDF and the reduced PDF of 
the bulk sample. From eqs. (1), (2), (23) and (24), the PDF and 
the reduced PDF are represented as 

8(r) = ^f(r) goo (r) + l-f(r), (25) 
Po 

G{r) = f(r)G QC (r)+47rr(p'o - Po)f(r). (26) 

The RDF, the PDF and the reduced PDF of the nanoparticle are 
obtained from those of the bulk sample modified by f(r) and 
the ratio of p' Q /po. The above g(r) and G(r) satisfy the normal- 
ization relations of g{r) — > 1 and G(r) — ► in the limit of 
r —* oo because /(r) — > for r — > oo. Then the eqs. (23), (25) 
and (26) are valid, and they can be applied to the PDF analysis 
on nanoparticles. 

In eq. (26), the second term is independent of the arrange- 
ment of the atoms included in the nanoparticle, and it depends 
only on the shape and the size of the nanoparticle (and the den- 
sities), while the first term depends on both the atomic arrange- 
ment in the nanoparticle, and the shape and the size of the 
nanoparticle. As a result, in the general case, we can expect 
that the peak structure due to the atomic arrangement given 
by the first term is on the "back ground" due to the shape of 
the nanoparticle given by the second term which can be eas- 
ily calculated by using f(r). From eq. (26), the "back ground" 
is larger for smaller po- Then the detailed experiment on the 
diluted nanoparticle can confirm the above expectations. 

Here, we discuss the difference between the PDF analysis 
and the small angle scattering. Generally, the shape and the size 
of the nanoparticles are determined by the small angle scatter- 
ing, and the scattering functions for the particles with various 
shapes have been calculated. However, since in the small angle 
scattering measurement, the wave vector is much smaller than 
the inverse of the atomic scale and the particle can be regarded 
as a continuum, the scattering profile depends on the shape and 
the size of the particle and is independent of the atomic arrange- 
ment. It indicates that the small angle scattering can not dis- 
tinguish the local disordered lattice from the ordered lattice, 



although, generally, they coexist in the nanoparticle. On the 
other hand, the PDF analysis can distinguish the local disor- 
dered lattice from the ordered lattice, because the functions of 
the PDF analysis depend on the atomic arrangement and the 
coherence. Then the PDF analysis can discuss the local lattice 
disorder in the nanoparticle by considering the finite size effect 
of the nanoparticle. Gilbert et al. (2004) have shown that the 
PDF profile of the zinc sulfide nanoparticle decreases with r 
more rapidly than the profile expected from the shape and the 
size of the particle. They pointed out that such reduction of 
the PDF profile is due to the local structural disorder driven 
by the strain in the particle and estimate the distance in which 
the structural coherence remains, from the reduction of the PDF 
profile. They speculate that such strain is caused by the irregular 
surface. 

Here, we can present a method to estimate the sizes of the 
domains with the ordered and the disordered lattice in nanopar- 
ticles, based on the present calculations. In the PDF analy- 
sis, the atomic density in the disordered part becomes nearly 
a continuum. For the continuum, Goo (r) is zero, which can be 
obtained from the relation /?oo(r) = 4Trr 2 p' Q , and eq. (24). Then 
the reduced PDF of continuum with nano-size is given by the 
smooth curve which corresponds with the second term of eq. 
(26). The reduced PDF of the part with the structural coherence 
shows sharp peaks corresponding to the atomic arrangement. 
Let's consider a spherical nanoparticle with a disordered sur- 
face, as a simple example. In this case, we consider two part 
: the sphere with a periodic atomic arrangement which has a 
diameter smaller than that of the particle, and the spherical shell 
of a continuum as shown in Fig. 10(b). The reduced PDF of the 
former are represented by modifying the functions of the bulk 
sample by using eq. (26). G(r) of the latter is given by the sec- 
ond term of eq. (26), obtained by using f(f) of the spherical 
shell with a corresponding size. We may calculate a sum of the 
two kinds of the reduced PDFs and fit it to the observed data 
by adjusting the thickness of the disordered surface, and as a 
result, we can estimate thickness of the disordered surface. By 
using our calculations, the information on the local inhomogen- 
ity in nanoparticles may be detected quantitatively. 



5. Conclusion 

The functions of the PDF analysis calculated for the various 
nanoparticles strongly depend on their shapes and sizes. The 
exact equations of the PDF analysis on the nanoparticle are pre- 
sented, by considering the correction factor f{r) and the ratio 
between the density in the particle and the total density aver- 
aged in the sample. The analysis which takes account of the 
finite size effect also enables us to estimate the sizes of the 
local disordered lattice and the ordered lattice included in the 
nanoparticle. 

The finite size effect on the PDF is remarkable in the large 
r-region where the accurate experimental data is hard to be 
obtained. The detailed analysis is enabled by a diffractometer 
with both of high resolution and high intensity, which can be 
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installed only in intense pulse neutron or intense synchrotron 
radiation source facilities. 
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